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Abstract 

Nanostructured materials exhibit low thermal eonduetivity beeause of the additional seat- 
tering due to phonon-boundary interaetions. As these interaetions are highly sensitive to the 
mean free path (MFP) of a given phonon mode, MFP distributions in nanostruetures ean be 
dramatieally distorted relative to bulk. Here we ealeulate the MFP distribution in periodie 
nanoporous Si for different temperatures, using the reeently developed MFP-dependent Boltz¬ 
mann Transport Equation. After analyzing the relative eontribution of eaeh phonon braneh 
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to thermal transport in nanoporous Si, we find that at room temperature optical phonons con¬ 
tribute 18% to heat transport, compared to 5% in bulk Si. Interestingly, we observe a steady 
thermal conductivity in the nanoporous materials over a temperature range 200 K < T < 300 K, 
which we attribute to the ballistic transport of acoustic phonons with long intrinsic MFP. These 
results, which are also consistent with a recent experimental study, shed light on the origin 
of the reduction of thermal conductivity in nanostructured materials, and could contribute to 
multiscale heat transport engineering, in which the bulk material and geometry are optimized 
concurrently. 


Introduction 

The quest for high-effieieney thermoeleetrie materials may be advaneed by using the ability of 
nanostruetures to suppress heat transport by several orders of magnitude with respeet to bulk with¬ 
out degrading eleetrieal transport signifieantly. ^ This phenomenon is based on the faet that phonon 
mean free paths (MFPs) are generally larger than eleetron MFPs; eonsequently, heat transport ex¬ 
hibits stronger size effeets. The extent of the suppression of phonon transport depends on the ratio 
between the intrinsie phonon MFP and the eharaeteristie length of the nanostrueture, L^. This 
ratio is known as the Knudsen number {Kn). When Lc is mueh smaller than MFP, i.e. for small 
Kn, phonon interaetions with boundaries are negligible. In this regime, heat transport reduetion is 
only due to geometrieal effeets, sueh as material removal in nanoporous materials, while phonon¬ 
boundary seattering is minimal. Therefore, heat transport is dominated by intrinsie seattering. On 
the other hand, for high Kn seattering is dominated by phonon-boundary interaetions. Within this 
regime, the phonon MFPs in the nanostrueture approaeh Lc and the phonons are eonsidered to 
travel ballistieally. The intermediate regime {i.e. Kn ^ 1) is often refered to as the quasi-ballistie 
regime. 

This analysis implieitly assumes single-phonon mode materials, but in most materials, there is 
a wide distribution of phonon MFPs, whieh in some eases span several orders of magnitude. For 
example, first-prineiples ealeulations for Si show that about half of the heat is earned by phonons 
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with MFPs larger than 1 jum. ^ Recent experimental measurements showing a reduction in thermal 
conductivity of Si membranes with microscale pores ^ provide support for these computational 
results. Together, they suggest that an accurate analysis of thermal transport in nanostructures 
should include the actual bulk MFP distribution. 

In bulk Si, the optical and acoustic phonons have very different MFP distributions. Optical 
phonons have relatively low MFPs because their dispersion curves are flatter than those of acoustic 
phonons, which by contrast have large MFPs. This effect has important consequences on thermal 
conductivity. First-principles calculations show that optical phonons contribute only 5% to the 
total thermal conductivity of Si, while the 95% is dominated by acoustic phonons.^ As a result, 
optical phonons are often neglected when calculating nanoscale heat transport in Si. However, in 
nanostructures, heat carried by optical phonons is slightly lowered while acoustic phonons can be 
strongly suppressed, making the two contributions comparable. 

A simple model for phonon-boundary scattering that neglects intrinsic scattering and assumes 
that phonons travel ballistically was devised by Casimir in 1938.^ Within the Casimir approach, 
the MFPs are considered to be the same as the characteristic length. By using the Casimir model, 
Tian et al^ concluded that optical phonons in Si nanowires contribute over 20% to the total thermal 
conductivity at room temperature. When dealing with complex boundaries, the Casimir approach 
fails for two reasons.^ First, it assumes that the characteristic length is known a priori, while in 
most materials with complex geometry this quantity is unknown. Secondly, a portion of the MFP 
distribution may lie in the diffusive or quasi-ballistic regime. 

In this work, we use the MFP-dependent Boltzmann Transport Equation (MFP-BTE)^ to cal¬ 
culate heat transport in nanoporous materials and provide the relative contribution of each phonon 
branch to the thermal thermal conductivity as a function of temperature. The use of the BTE 
enables treatment of complex geometries with a good level of predictive power. We focus on 
nanoporous Si (np-Si) with aligned pores with square cross section and a periodicity of 10 nm, and 
we consider the temperature range 100 K — 300 K. We show that at room temperature the thermal 
conductivity in np-Si is suppressed by more than one order of magnitude with respect to bulk Si, 
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with longitudinal optical (LO) phonons contributing nearly 20% to the total heat transport. This 
result is in agreement with the previous qualitative discussion. Further, we find that the thermal 
conductivity of np-Si exhibits a plateau over the temperature range 200 K — 300 K. We demon¬ 
strate that this arises from two effects: first, as most of the acoustic phonons travel ballistically 
because of their large Kn, their MFPs in np-Si are constrained by the characteristic length of the 
material. Second, in this temperature range, the heat carried by optical phonons changes weakly 
with temperature in bulk Si itself, inducing similar behavior to np-Si. By revealing the microscopic 
mechanisms leading to the reduction in heat conduction, our findings may enable new approaches 
for engineering high-efficiency thermoelectric devices. 


Results and discussion 


To compute the reduction of heat transport in nanostructures, we employ the concept of the “sup¬ 
pression function,” 5(A), which defines the departure from diffusive transport in terms of the MFP 
distribution function,^ 


5(A) 




( 1 ) 


where ^^“^'^(A) and .firp™°(A) are the bulk MFP distributions for branch p in bulk Si and np-Si, 
respectively. Within the relaxation time approximation, the effective thermal conductivity for each 
phonon branch can be written as 


^nano ^ / A';(««"(A) JA = / K^f^{A)S{K)dA. (2) 

^ ./o Jo ‘ 

The total thermal conductivity is then given by In the case of purely diffusive 

transport, the suppression function is MFP-independent and Eq. (2) leads to the diffusive thermal 
conductivity where g is a function that depends only on the material geometry and 

j^buik jg thermal conductivity. The bulk MFP distribution at different temperatures can 

be obtained either experimentally through MFP reconstruction techniques * or computationally. In 
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this work we adopt a first-principles approach based on Density Functional Theory (DFT) and 
the linearized BTE.^’^ The bulk MFP distribution is computed via where 

05^“"^ (A) is cumulative thermal conductivity. We recall that the cumulative thermal conductivity is 
the thermal conductivity of phonons whose MFPs are below a given A. The intrinsic MFP for 
a given phonon mode is defined as A = |v|t'’“*'^. We note that does not include boundary 

scattering. For this reason, in rest of this study, we will refer to such a result as bulk-BTE. 

The suppression function can be obtained in different ways, depending on the system and the 
required accuracy. In some cases, such in nanowires and thin films, 5'(A) can be obtained analyti¬ 
cally within a reasonable level of accuracy. However, most of the analytical derivations are based 
on the “gray” approximation, which assumes phonon dispersions described by a single group ve¬ 
locity. Eurthermore, formulae for the suppression function are limited to simple geometries. In 
this work we therefore employ a recently developed formulation of the BTE that requires only 
the bulk MEP distribution This method, MEP-BTE, has the same accuracy as the the 

commonly used frequency dependent approach (ED-BTE), provided that we consider small ap¬ 
plied temperature gradients, AT/L, where AT is the applied difference of temperature and L is 
the distance between the hot and cold contact. The key equation of the MEP-BTE consists of the 
integro-differential equation 

As-Vf(A)-ff(A) = r£ / ^l^^<f(A')>JA', (3) 

where T (A) represents the normalized temperature associated with phonons having MEP A, given 
by T (r, A) = In Eq. (3), s is the phonon propagation direction and 7 a material property 

y^bulk ^ 

given by 7 = '£pl^ -^dA = 2.2739 x 10^^^ KW^^m^^ for Si. The notation < .r > stands 
for an angular average. The right-hand-side of Eq. (3) is equal to where is the effective 

lattice temperature,^ which does not depends on A explicitly, and provides an average of the local 
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energy of phonons. Once Eq. (3) is solved, the suppression function can be computed via 


S{A) = ^J^<fs-ns>dx^, (4) 

where E is the surface of the hot contact having normal and area A. 

Our simulation domain consists of a square unit cell, which we choose to have a size of 
L = 10 nm. The unit cell, which contains one square pore, is subjected to a difference of temper¬ 
ature AT = 1 K about Tq. Periodic boundary conditions are applied to both the longitudinal and 
transverse direction of heat flux, nf, which is enforced by applying a difference of temperature AT 
along nf, i. e. 

T (A, s, r) - f (A, s, r + P) = (n ■ nf) AT, (5) 

where r runs along the faces of the unit cell, n is the normal to the boundary pointing outside 
the domain and P is the periodicity vector. It is straightforward to show that along the direction 
perpendicular to the heat flux, no difference of temperature is imposed. We assume that we have 
an infinite material along the directions orthogonal to the pore plane. Moreover, we assume that 
phonons are scattered diffusively when interacting with the pores, i.e., the distribution of phonons 
moving away from a pore’s surface is uniformly distributed in angular space. We note that this 
assumption must be used with caution, especially at very low temperatures, as some phonons can 
be reflected specularly depending on the roughness of the boundary. However, in this work we 
assume the surfaces have significant roughness such that specularity effects can be neglected in the 
considered temperature range. Details on the surface specularity effects on thermal transport can 
be found in Ref. 

We first validate the code with the case of zero porosity, i.e., bulk Si. Figure Figure 1 illustrates 
the good agreement between experimental data, bulk BTE predictions, and MFP-BTE calcula¬ 
tions of the thermal conductivity of bulk Si. The bulk case recovers purely diffusive transport, as 
there are no boundaries suppressing phonons. 

Figure Figure 1 also shows the thermal conductivity for np-Si with porosities (j) = 0.05 and 
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0 = 0.25. At r = 300 K, we observe a reduction of more than one order of magnitude with 
respect to the bulk, in accordance with our previous studies.We now analyze the relative con- 
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Figure 1: The thermal conductivity for different temperatures of bulk Si and np-Si. For bulk Si, the 
experimental data, bulk-BTE and MFP-BTE calculations are shown and agree with each other. 
For np-Si, we show the cases of porosity (j) = 0.05 and (j) = 0.25, computed by the MFP-BTE. In 
the inset, we show the normalized temperature map of the np-Si with porosity (j) = 0.25. 


tribution of each phonon branch to the total thermal conductivity. In Eig. Eigure 2-a, the normalized 
cumulative thermal conductivity for bulk Si at T = 300 K, computed by the the bulk BTE, is shown. 
As expected, acoustic phonons contribute most to the thermal conductivity, while optical phonon 
contributions are small. In particular, the two transverse acoustic (TA) branches and longitudinal 
acoustic (EA) branch contribute approximately one third each to the total thermal conductivity. 
The figure also shows that the longitudinal optical (EO) phonons contribute 5% to thermal trans¬ 
port, while the transverse optical (TO) phonons have a negligible contribution. However, the EO 
phonons start to contribute significantly in np-Si, reaching 16% of the total thermal conductivity 
in the case of 0 = 0.05 (Eig. Eigure 2-b). The TO contribution remains negligible. 

The roughly four-fold increase in the relative contribution of EO phonons can be better under¬ 
stood by analyzing the MEP distributions in relation to L^.. According to Ref., in porous materials 
Lc can be defined as the pore-pore distance in the direction orthogonal to thermal flux. The pore- 
pore distance in an array of square aligned pores is related to the porosity via Lc — L{ \ — a/^). 
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which leads to the values 5 nm and 7.76 nm for ^ = 0.25 and ^ = 0.05, respectively. The charac¬ 
teristic length dictates the transport regime of phonons with a given MFP. Figure Figure 2-a shows 
that the maximum MFP of LO phonons contributing to the thermal conductivity is around 20 nm, 
while acoustic phonons have MFPs up to 10 /im. As a result, optical phonons, which generally 
have MFPs similar to Lc, are less suppressed than acoustic phonons. For (j) = 0.25, the characteris¬ 
tic length is even smaller and, consequently, the relative LO phonon contribution increases (up to 
18%, as shown in Fig. Figure 2-c). 

When Lc is larger (e.g., 100 nm), the effect of the nanostructure on optical phonons becomes 
negligible, but most acoustic phonons are still suppressed. In this case, it is possible to have a 
“reversal effect,” in which optical phonons are the main contribution to the thermal conductivity. 

For macroscopic samples, e.g. Lc > 100 /im, the thermal conductivity approaches the value 
predicted by the Fourier model and the MFP distributions are restored to the bulk ones times the 
geometric factor, g, that depends only on the geometry. For aligned porous materials, the geometry 
factor can be well approximated by g = This approximation was validated against finite- 

element modeling of diffusive heat conduction. 

This finding has important consequences on the guidelines for optimizing nanostructured ther¬ 
moelectric materials. Typically, the geometry of the nanostructure and the bulk thermoelectric 
materials are optimized separately. Here we suggest that both macro and nanoscale have to be 
considered concurrently. The following example helps clarify this point. Let us assume that we 
have two “gray” materials, A and B, with average MFPs and As, respectively. We further as¬ 
sume that the thermal conductivity of material B is larger than that of material A. We consider a 
nanostructure with Aa « Lc << Ag. Material B will undergo strong phonon suppression whereas 
heat transport in material A will still be in the diffusive regime. It is clear therefore that, with a suf¬ 
ficiently large Ag, material B exhibits lower thermal conductivity than that of material A, resulting 
more appealing for thermoelectrics. Similar conclusions can be drawn for non-gray materials. 

We now investigate the temperature dependence of thermal conductivity in the range 200 K < T < 300 K. 
All the results shown below refer to the case with (j) = 0.25. Similar conclusions can be drawn for 
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Figure 2: Cumulative thermal conductivity at T = 300^ for a) bulk Si, b) np-Si with porosity 
(j) = 0.25 and c) np-Si with porosity (j) = 0.05. All the values are normalized to the total thermal 
conductivity. The largest MFP contributing to heat transport for a given branch can been seen from 
the point where the relative cumulative thermal conductivity becomes flat. 


the case with 0 = 0.05. As shown in Fig. Figure 1, the thermal conductivity of np-Si exhibits 
little change in this temperature range, whereas it decreases as 1/T due to Umklapp scattering in 
bulk Si.^’^ This behavior arises from the very large Kn of acoustic phonons, which therefore travel 
ballistically. Within this regime, the phonon suppression function goes as ^(A) Lc/A, and the 
contribution of acoustic phonons in np-Si becomes independent of the bulk MFP. According to 
Fig. Figure 3, for temperatures as low as 200 K, the contribution to the total thermal conductivity 
from large MFP acoustic phonons becomes even greater. Optical phonons, on the other hand, have 
MFPs close to Lc, and, in principle, their temperature dependence in bulk Si would affect their 
MFPs in np-Si. However, according to Fig. Figure 3-a, for temperatures 200 K < T < 300 K, 
heat carried by LO phonons in bulk Si does not change significantly with temperature because the 
increase in heat capacity is compensated by the decrease in scattering time.^ Consequently, heat 
carried by LO phonons in np-Si does not change with temperature. Heat carried by TO phonons 
is negligible in both bulk Si and np-Si. These combined effects lead to the observed insensitivity 
of thermal conductivity to temperature in np-Si. Finally, as the temperature decreases in the range 
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Figure 3: MFP distributions of bulk Si at a) T = 200 K and b) T = 300 K. In both panels, the 
shaded area shows the region of MFPs that are elose to the eharaeteristic size of the nanostrueture, 
Lc = 5 nm. At low temperatures, the contribution to the thermal conductivity from long-MFP 
acoustic phonons rises. However, the ballistic regime constrains these MFPs to be equal to Lc. 


100 K < T < 200 K because of the decrease of the heat capacity, heat carried by LO phonons in 
bulk Si starts to decrease, as shown in Fig. Figure 4-a. As a consequence, their relative contribu¬ 
tions to the thermal conductivity in np-Si decreases, as well. 

Our prediction of the temperature-dependence of thermal conductivity has been confirmed ex¬ 
perimentally in a recent study, where heat transport has been measured in holey Si with periodicity 
of about 20 nm. We point out that our model, being based on BTE, calculates phonon dynam¬ 
ics within the particle picture, while the wave picture is only retained in the calculations of the 
bulk dispersion curves. As a consequence, our study suggests that in such systems incoherent 
scattering dominates thermal transport. The importance of coherence effects has been assessed in 
another recent study, where Monte Carlo simulations were used to compute thermal transport in 
np-Si membranes. Their conclusion was, however, that incoherent effects can fully explain the 
remarkably low thermal conductivity in silicon nanomeshes. 
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Figure 4: Contribution of each phonon branch to the total thermal conductivity for a) bulk Si and 
b) np-Si. In the range 200 K < T < 300 K, the contributions to the thermal conductivity from each 
phonon branch in np-Si do not change significantly with temperature. 


Conclusions 

Using the MFP-BTE, we calculate the temperature dependence of thermal conductivity in np-Si. 
We quantify the contribution of optical phonons to thermal conductivity in np-Si with periodicity 
10 nm, which at room temperature amounts to 18%. We also predict constant thermal conductivity 
over the range 200 K < T < 300 K, which has been recently observed experimentally. Our findings 
help further the understanding and manipulation of heat transport at the nanoscale for low thermal- 
conductivity applications such as thermoelectrics. We have also showed that the effectiveness 
of nanostructuring in reducing thermal transport does not depend on the bulk thermal conductivity 
but rather on the bulk MFP distributions of phonon branches. Consequently, our approach suggests 
that the search for better nanostructured thermoelectric materials has to involve the shape of the 
bulk cumulative thermal conductivity in relation with the material’s geometry. In other words, the 
material optimization has to be done at both macro and nanoscale concurrently. 


11 



















Methods 


We define the eumulative thermal eonduetivity in the bulk material as 

^buiM = ( 6 ) 

where Cp(q) is the heat eapacity, Vp,x(q) is the group veloeity along the ;c-direetion, Tp“^'^(q) is 
the three-phonon seattering time, and 0 is the Heaviside funetion. We reeall that the eumulative 
thermal eonduetivity is the thermal eonduetivity of phonons whose MFPs are below a given A. 
All the quantities appearing in Eq. (6) are taken from Ref.^ and are not reported here for the 
sake of simplicity. The phonon dispersion curves and scattering times are obtained by means of 
harmonic and anharmonic force constants, which are extracted from DFT. The system’s relaxation 
times are computed by using a uniform reciprocal space grid of 24 x 24 x 24 points, harmonic 
force constants up to 5'* neighbors and cubic force constants up to first neighbors. We use the 
Local Density Approximation (EDA) from Perdew and Zunger^^ with a energy cutoff of 40 Ryd. 
The MFP-BTE is solved for a set of 30 MFPs, uniformly spaced on logarithmic scale from about 
0.1 nm to 100 jum. The MFP-BTE is discretized both in angle and space. The spatial discretization 
is achieved by means of the finite volume (FV) approach whereas the solid angle is discretized by 
means of the discrete ordinate method. 
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